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Charged Current Single Pion Cross Section Measurement at MiniBooNE 

M.O. Wascko a , for the MiniBooNE Collaboration 

a Department of Physics and Astronomy, 
Louisiana State University, Baton Rouge, LA 70803 

We present MiniBooNE's preliminary ^ m CC1tt + cross section measurement, calculated using the ratio of CCl7r + 
to CCQE events. We find the inclusive CCl7r + measurement to be below the nuance pQ and NEUGEN [5] 
expectations. 
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1. Introduction 

Charged current single pion (CC1vt + ) produc- 
tion has been studied since the advent of high en- 
ergy neutrino beams, but the cross section around 
1 GeV energy is still not well understood. Also, 
many of the data that do exist come from hy- 
drogen and deuterium targets, so there is still a 
significant need to study nuclear effects in this 
process. Figure ^ |3J shows a comparison of the 
CCl7r + measurements from the Argonne U and 
Brookhaven |S] bubble chamber experiments, as 
well as the nuance [I] Monte Carlo prediction for 
the CC1vt + cross section. 

The MiniBooNE CCf7r + event sample is a 
semi-inclusive sample because the Cherenkov 
calorimeter detector is not able to resolve the fi- 
nal state recoil nucleons, and hence cannot dis- 
tinguish the exclusive channels. In this work, we 
use the label CCi7r + to indicate either of the res- 
onant reactions, v^p — ► /i~A ++ — * /i~p7r + or 
Vfjji — > (i~ A + — > /i~mr + , as well as coherent pion 
production, v^A — > /i - Att + . A significant frac- 
tion of the CCl7r + events at these energies are 
expected to arise from coherent production 6 . 
Given the importance of CCl7r + events as a back- 
ground for — ► v x disappearance searches |7], 
the recent K2K SciBar limit on CC coherent pion 
production at 1.3 GeV [Sj, and the importance of 
coherent neutral current tt° (NC7t°) production 
as a background source for — » v e oscillation 
searches 0, there is significant interest in mea- 
suring the CCl7r + cross section on carbon. 

MiniBooNE |U| is a neutrino oscillation experi- 
ment at Fermilab designed to confirm or rule out 
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Figure 1. Previous measurements of the total 
cross section per nucleon of the process v^p — > 
pTpTi + at low neutrino energy. 



the hypothesis that the LSND v e excess ^0] is due 
to — ► v e oscillations. A general description 
of the experiment can be found elsewhere 
CCl7r + events are expected to comprise ~25% of 
the total MiniBooNE neutrino event rate, mak- 
ing these the second most probable interactions 
after charged current quasi-elastic (CCQE). For 
this reason they are interesting and useful for 
MiniBooNE. Prior measurements of these cross 
sections suffered from poor statistics ^21- Sta- 
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tistical precision will not be a problem for Mini- 
BooNE; the MiniBooNE Monte Carlo, which uses 
the nuance neutrino generator ,T|, predicts that 
we should have almost 58,000 i^CCl7r + events af- 
ter cuts in the MiniBooNE detector with the full 
data set of 5.5xl0 20 protons on target (POT). 

1.1. Cross Section Ratio Overview 

In this analysis, we normalize the observed rate 
of CC17T+ events to that of CCQE events, and 
equate that to the ratio of CC17T+ to CCQE 
cross sections. In doing the analysis this way, 
we use the fact that the same neutrino flux gen- 
erates both event samples. Using the ratio al- 
lows us to neglect the uncertainties in the neu- 
trino flux prediction. Moreover, if disappear- 
ance were present in the data, the predicted num- 
ber of events in the Monte Carlo would be incor- 
rect because of the depletion of the flux due 
to Vfj, — y v x oscillations. By normalizing to the 
CCQE data we avoid this issue. 

We write the CC17T+/CCQE cross section ratio 
as: 



o~ccqe{E v ) 



(1) 



where N® ata is the true number of events of type 
a in the data. The true neutrino energy is de- 
noted by E v . If we assume that the CCQE cross 
section is well simulated by the Monte Carlo, then 
we may rewrite Equation^ 
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and thus obtain the CCl7r + cross section mea- 
surement as a function of neutrino energy. 

The raw event ratio is measured in the Monte 
Carlo as: 

(E? EC ) 
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where we have used E^ EC to denote the recon- 
structed energy of events from either process. To 
equate this to the cross section ratio, we must 
account for energy smearing, cut efficiencies, and 
the presence of background events in each data 
sample. We use the Monte Carlo to estimate each 
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Figure 2. Distribution of events versus re- 
constructed neutrino energy, for events passing 
CCl7r + cuts. Black points show data with statisti- 
cal errors. Monte Carlo results are shown in yel- 
low, with error bars showing some important con- 
tributions to the systematic uncertainties. Data 
and Monte Carlo are normalized to unit area. 



of these factors, and apply the derived corrections 
to the data samples before equating the event ra- 
tio to the cross section ratio. 

2. CC17T+ Event Analysis 

MiniBooNE's CClvr+ event selection requires 
the simple yet robust cut of two Michel electrons 
following the neutrino interaction Approx- 
imately 40% of pions emitted at these energies 
stop in the detector oil. These decay to muons 
with lifetime r = 2.6 x 10 -8 ns, which then 
decay to Michel electrons. The muons emit- 
ted from the neutrino interaction also come to 
rest, and 92% of these decay to Michel electrons. 

Applying this requirement to 3.3xl0 20 POT 
of MiniBooNE data yields over 44,000 CC1tt + 
candidate events, making this data set larger by 
a factor of five than all previous CCl7r + bub- 
ble chamber data published to date. The Monte 
Carlo predictions indicate that the event selection 
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Table 1 

Event composition of CCl7r + sample. 



Reaction Type 


Percentage 


resonant CC17T+ 


75.8% 


coherent CCl7r + 


9.2% 


CC QE 


4.1% 


multi-pion 


6.1% 


DIS 


2.6% 


CCvr 


1.5% 


other 


0.7% 



cuts are ~30% efficient for CCl7r + events within 
the fiducial radius of 500 cm, with a purity of 
85%. Tabled shows the fractions of signal and 
background events passing the CCl7r + event se- 
lection. The background events come from either 
events with multiple pions that lose one or more 
of them within the nucleus, or CCQE events that 
acquire a 7r + though hadronic interactions within 
the nucleus, so that the events contain a single 
fx~ and 7r + in the final state. 

The CCl7r + events are currently reconstructed 
under simple assumptions: we apply a single ring 
fitter to the PMT hits to find the position, di- 
rection, and photon flux from the Cherenkov ring 
produced by the \x~ in the event. This yields the 
muon energy using only the Cherenkov light in 
the reconstructed ring; this is done to avoid in- 
cluding light generated by the ir + in calculating 
the muon energy. We then use the fitted energy 
and direction of the muon to reconstruct the neu- 
trino energy by assuming the CCl7r + reaction is 
a simple two body collision, e.g. v^p — > /i~A ++ . 
Then we write the neutrino energy as: 
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where m p is the proton mass, is the muon 
mass, to a is the mass of the recoil resonance 
(which we take to be exactly 1232 MeV), and 
and are the reconstructed energy and angle of 
the outgoing muon. 

Figure |21 shows the fraction of CC17T+ events 
as a function of reconstructed neutrino energy, 
for both data and Monte Carlo. The data are 
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Figure 3. Distribution of muon lifetimes from 
charged current single pion events. Only data are 
shown, with only statistical errors. 



shown as black points, with error bars represent- 
ing the statistical uncertainties. The Monte Carlo 
results are shown as yellow points, with yellow er- 
ror bands representing the systematic uncertain- 
ties. Monte Carlo studies indicate that the neu- 
trino energy resolution is about 20%. Much of 
the resolution comes from the assumption that 
the recoil resonance has mass equal to 1232 MeV; 
however MiniBooNE is not able to reconstruct 
the invariant mass of the A because the emitted 
nucleon does not leave a Cherenkov ring. The 
CCl7r + Monte Carlo predictions are based on the 
Rein and Sehgal model for resonant and coherent 
pion production 6 14j, and the nuance nuclear 
model pQ. 

The Michel electrons from the CCl7r + candi- 
date events are used to verify the composition 
of the data set. The fi~ are captured by car- 
bon nuclei with a probability of 8%, changing the 
lifetime to 2026. 3±1. 5 ns ^H] from the expected 
2197.03±0.04 ns EE]. The distance from each re- 
constructed Michel to the end of the fi~ track 
is calculated, the Michels are sorted into "close" 
and "far" samples based on the distance to the 
end of the /i~ track. The lifetimes for the two 
samples are shown in Fig. EI The observed muon 
lifetimes for the close and far Michel samples are 
2057±14 ns and 2218±15 ns, respectively, indi- 
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eating that they are indeed mostly from /i - and 
fi + decay, respectively. The same studies per- 
formed on the Monte Carlo events yield lifetimes 
of 2043±15 ns and 2203±15 ns, and show that 
the close and far samples are ~80% pure \i~ and 
fi + , repectively. 

3. CC QE Event Selection 

The CCQE event reconstruction is described 
in greater detail elsewhere |17j . To summarize, 
the event selection uses a 10 variable Fisher dis- 
criminant designed to find a single Cherenkov 
ring from a muon. The neutrino energy is recon- 
structed using simple two-body kinematics and 
corrected for the effects of Fermi momentum of 
the target nucleon. Figure^lshows the fraction of 
CCQE events as a function of reconstructed neu- 
trino energy, for both data and Monte Carlo. The 
data are shown as black points, with error bars 
representing the statistical uncertainties. The 
Monte Carlo results are shown as yellow points, 
with yellow error bands representing the system- 
atic uncertainties. 

Applying the CCQE event selection to 
2.3 xlO 20 POT of MiniBooNE data yields over 
60k CCQE events, with 86% purity. Table 
shows the fractions of signal and background 
events passing the CCQE event selection. The 
background events come mostly from single pion 
events in which the pion is absorbed within the 
nucleus, so that the events contain a single pT 
and no 7r in the final state. 

Monte Carlo studies indicate that the neutrino 
energy resolution is about 10%. This is bet- 
ter than the CCl7r + energy resolution because 
the mass of the recoil particle (a proton) is well 
known, and the muon Cherenkov ring is cleaner 
because there is no n + in the final state. 

4. CC17T+ Cross Section 

To use Equation [21 to calculate the cross sec- 
tion, we first calculate the ratio of CCl7r + /CCQE 
events. The numerator of the ratio is the event 
sample shown in Figure after it was corrected 
for the background fraction, energy smearing and 
cut efficiencies. Similarly, the denominator is 



Table 2 

Event composition of CCQE sample. 



Reaction Type 


Percentage 


CCQE 


86.0% 


CCl7r + : resonant 


8.9% 


CCl7r + : coherent 


1.4% 


NC17T+ 


1.7% 


CC17T 


1.1% 


CC^° 


1.5% 


other 


0.9% 



shown in Figure H The ratio of CC1tt+/CCQE 
events as a function of neutrino energy is shown 
in Figure [5] Note that this neutrino energy has 
been corrected for energy smearing. The error 
bars show the systematic uncertainties added in 
quadrature with the statistical uncertainties. 

We restrict the analysis to the reconstructed 
energy region 0.5 GeV-1.4 GeV. This is because 
the CC17T+ sample has insufficient statistics be- 
low ~0.5 GeV, due to the energy threshold for 
A production, and the CCQE sample has insuf- 
ficient statistics above ~1.4 GeV, due to the effi- 
ciency of the event selection cuts for high muon 
energies. 

To calculate the CC17T+ cross section, we mul- 
tiply the ratio by the Monte Carlo CCQE pre- 
diction. The Monte Carlo prediction is based on 
the nuance neutrino generator, which assumes a 
Llewellyn-Smith quasi-elastic cross section |18| . 
and a Smith-Moniz Fermi gas model for nuclear 
interactions |19|. The nuance program also sim- 
ulates final state interactions, which mainly in- 
clude the interaction probabilities of the final 
state particles while they are still inside the nu- 
cleus 0. The MiniBooNE Monte Carlo uses the 
BB A03 non-dipole vector form factors [20j , and a 
dipole axial form factor with Ma=1-03. 

Figure shows MiniBooNE's CCl7r + cross sec- 
tion measurement. Again, the data are shown as 
black points, with statistical error bars. The yel- 
low error bands represent the present estimate 
of the systematic uncertainties. Also shown are 
the nuance and NEUGEN 2 predictions for the 
CCl7r + cross section. Although the CCl7r + mea- 
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Figure 4. Distribution of events versus re- 
constructed neutrino energy for events passing 
CCQE cuts. Black points show data with statisti- 
cal errors. Monte Carlo results are shown in yel- 
low, with error bars showing some important con- 
tributions to the systematic uncertainties. Data 
and Monte Carlo are normalized to unit area. 



Figure 5. Ratio of CCln + /CCQE events vs. 
neutrino energy, after corrections for background 
events, energy smearing, and cut efficiencies. 
The black error bars show the statistical and sys- 
tematic uncertainties added in quadrature. 



surement is calculated from the ratio with the 
nuance CCQE cross section, the comparison to 
NEUGEN is appropriate because its CCQE pre- 
diction is identical to the nuance prediction. 

Note that the MiniBooNE CC1tt + measure- 
ment lies ^25% below the nuance and NEUGEN 
predictions. One plausible explanation for this 
apparent disparity is offered by Figure ^ the 
AM. Jj and BNL[5J measurements of the v^p — > 
/j,~Ptt + cross section. The BNL measurement lies 
~40% above the ANL, and nuance splits the dif- 
ference. MiniBooNE's measurement lying 25% 
below the nuance prediction suggests that it is 
more consistent with the ANL observation. 

5. Systematic Uncertainties 

The sources of systematic error considered here 
are summarized in Table and can be gener- 
ally categorized as neutrino cross section, optical 
model, and energy scale uncertainties. By optical 



model, we mean the properties of light genera- 
tion and transmission in the detector mineral oil. 
The systematic uncertainties are assessed on the 
CC17T+/CCQE ratio, so there is some cancella- 
tion of errors. 

We include cross section uncertainties in a cross 
section measurement because we must assess the 
error on the predicted number of events which 
pass the CCQE selection cuts. The CCQE back- 
ground after cuts is dominantly from CCl7r + 
events, and so any normalization uncertainty on 
the predicted CCl7r + (and CCQE) event rates 
must be included. Second, efficiency corrections 
are applied to the CCl7r + data which are derived 
from the Monte Carlo, and so any energy depen- 
dent uncertainties must be included. 

The cross section uncertainties are all derived 
from external data [21] ■ The optical model un- 
certainties are derived from "table-top" measure- 
ments of the oil and in situ calibration data sam- 
ples |22] ■ The correlations between the sources of 
uncertainty are not yet fully determined, so we 
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Figure 6. Semi-inclusive CCl7r + cross sec- 
tion vs. neutrino energy, extracted assuming 
the nuance CCQE cross section prediction with 
yfj^ — 1.03 GeV. Statistical uncertainties are 
shown by the black error bars, and the impor- 
tant contributions to the systematic uncertainties 
are shown by the yellow error bands. Also shown 
are the nuance and NEUGEN predictions for the 
semi-inclusive CCl7r + cross section vs. recon- 
structed neutrino energy. 



present the analysis assuming they are all uncor- 
rected, adding all uncertainties in quadrature. In 
fact, we expect many of the sources of uncertainty 
to be anti-correlated. 



6. Conclusions 

We have presented the first measurement of the 
semi-inclusive h>^CClTr + cross section ratio on a 
nuclear target near 1 GeV. From this, we ex- 
tract a CCl7r + cross section measurement which 
is lower than the predictions of the nuance and 
NEUGEN Monte Carlos, but appears consistent 
with the ANL observation of resonant single pion 
production. 



Table 3 

Sources of uncertainty considered for the CCl7r + 
cross section analysis, summed over the entire 
analysis region from 0.5 to 1.4 GeV. 

source effect on a 



Optical Model 


20% 


Cross Sections 


15% 


Energy Scale 


10% 


Statistics 


6% 
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